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Abstract — An anisotropy signal for the arrival directions of 
ultra-high energy cosmic rays (UHECR) of more than 99% 
confidence level was established using data collected by the Pierre 
Auger Observatory. Cosmic rays with energy above ~ 6 x 10^^ 
eV show a correlation with the positions of extragalactic nearby 
active galactic nuclei (AGN), being maximum for sources at 
less than ~100 Mpc and angular separation of a few degrees. 
The evolution of the correlation signal with the energy shows 
that the departure from anisotropy coincides with the flux 
suppression observed in the spectrum, being therefore consistent 
with the hypothesis that the correlated events have their origin in 
extragalactic sources close enough to avoid significant interaction 
with the cosmic microwave background (the Greisen-Zatsepin- 
Kuz'min effect). Even though the observed signal cannot un- 
ambiguously identify AGNs as the production sites of UHECRs, 
the potential sources have to be distributed in a similar way. A 
number of additional statistical tests were performed in order to 
further understand the nature of the correlation signal. 

I. Introduction 

The origin of the highest energy cosmic rays (above 10^^ 
eV), as far as their mass composition and production mech- 
anisms are concerned, is still undetermined. However, the 
recent Auger results, which include limits on the diffuse flux 
of neutrinos [1], bounds on the flux of photons above 10^^ 
eV obtained with the Surface Detector (SD) and Fluorescence 
Detector (FD) together [2] and with the SD only [3], have put 
stringent bounds on top-down models, favoring, therefore, the 
hypothesis of production by acceleration process in powerful 
astrophysical sources [4]. 

Anisotropy studies play an important role in the identifica- 
tion of the sources and can even provide valuable information 
on the chemical composition. Detectors like the large Auger 
surface array are sensitive to showers with energies around 
above 10^^ eV, but the detector reaches 100% efficiency 
only above 3 x 10^^ eV. Given the steep falling spectrum, 
even taking into account the rapidly falling efficiency of the 
detector below this threshold, the Auger Surface Detector 
sample is dominated by below- saturation showers. In this 
energy region, despite the higher statistics, the myriad of 
systematic effects associated to the detector exposure makes 
large scale anisotropy analysis a very challenging task [5]. 
Weather induced effects, like seasonal and diurnal trigger rate 
modulations [6], geometric effects associated to the particular 
arrangement of the SD array and the influence of the local 
geomagnetic field on the air shower development are some of 
the examples of systematic effects which one has to deal with 



in order to extract reliable information on anisotropics over 
extended regions of the sky. For energies above 10^^ eV, the 
Pierre Auger Collaboration has already published some results 
on one-dimensional analyzes, that is, right-ascension first 
harmonic analysis, Fourier development in modified sidereal 
time [7] and the East- West method [8]. The galactic centre 
is a natural candidate for searches of anisotropy, since there 
are a number of evidences that it harbors, as essentially every 
galaxy in the universe, a massive black hole, which would be 
associated to the strong radio emissions from Sagittarius A*, 
and the 7-rays detected by H.E.S.S. [9]. However, above lO^*" 
eV, the Auger results in the galactic centre region show no 
signs of anisotropy [10], despite claims in the past of large 
excesses in this region of the sky [11], [12]. 

At the very edge of the cosmic ray spectrum measured so 
far, the systematic effects are no longer so dominant, however, 
now the extremely low statistics is the challenge to be faced. 
On the other hand, it has been long since people realized 
that at such extreme energies, charged particles magnetic 
rigidity, both in galactic and extra-galactic fields, should be 
large enough to avoid a complete isotropisation of their 
arrival directions. In fact, the Auger Collaboration has recently 
demonstrated that its highest energy events, that is, above 60 
EeV (1 EeV = 10^^ eV), are anisotropic [13], [14]. The result 
has been established by means of a correlation found between 
the arrival directions of the showers and a catalog of AGNs. 
The result has demonstrated the feasibility of charged particle 
astronomy. This conference paper contains a discussion of the 
results presented in the two papers mentioned above and some 
remarks on their implications to the cosmic ray physics as well 
as on issues they have raised after publication. 

II. Current Array Status and the Data Set Used 

The Pierre Auger Observatory has now reached completion 
with all of its 1600 surface detector stations deployed on 
ground and the FD with all 24 telescopes taking data together 
since February 2007. The surface detector has been taking 
data stably since January 2004. The data presented in this 
conference contribution corresponds to showers detected by 
the SD in the period 1 January 2004 up to 31 August 2007, 
amounting for a total integrated exposure of 9000 km^ sr yr. 
In this period, the surface detector has continuously grown in 
size from only 154 to 1388 stations. 
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Fig. 2. (a) Distribution of the space angle between reconstructions of the same shower performed with independent sets of tanks in a sub-array of doublets. 
Defining the AR as the angle containing 68% of all reconstructed showers coming from a point source, AR=1.5cr. (b) Distribution of the space angle between 
the reconstruction done using only the surface detector and the surface and fluorescence detectors together (the hybrids). Notice that the intrinsic hybrid 
resolution contributes to the total dispersion observed in this plot. 




Fig. 1. Surface detector array status as of 27 October 2008, showing a 
total of 1638 stations deployed and taking data. The field of view of the four 
fluorescence eyes (with 6 telescopes each) are represented by the blue lines 
(each one is 20 km long). Some infills, tanks deployed at half the distance 
of the regular array from each other in order to lower the observatory energy 
threshold can also be seen close to one of the eyes. 



In order to avoid border effects, which tend to be important 
during the detector growing process, when the array has a non 
negHgible number of holes, we consider only events where the 
station with the highest signal is surrounded by at least 5 other 
active tanks, that is, stations actually taking data when the 
shower reaches the ground. To further constrain the impact of 
dead stations in the event reconstruction, the estimated shower 
core position is required to fall inside a triangle of active 
stations. 

Obviously, for anisotropy studies, the accuracy with which 
one reconstructs the shower direction plays a special role. 



The reconstruction of the event direction is done by fitting 
the arrival times predicted by a certain shower front model 
(either a plane or a curved form) propagating at the speed 
of light to the measured arrival times in the tanks triggered 
by a shower. The accuracy achieved in such a reconstruction 
is limited essentially by the clock resolution of the station 
and the intrinsic fluctuations of the arrival time of the first 
particle [15]. The Auger SD angular resolution (AR) above 
10^^ eV, is estimated to be better than 1°, as defined by the 
angular window which would contain 68% of all reconstructed 
showers produced by primaries coming from a point source 



[16]. Figures |2(a)| and |2(b)| show two different estimations 
of the AR in Auger, using a sub-array of doublet^ through 
which is possible to reconstruct the same shower with two 
independent sets of tanks and using the hybrid events, by 
comparing the SD and SD+FD reconstructions. 

By profiting from the unique hybrid nature of the Auger 
detector, the so called hybrid events, the ones detected simul- 
taneously by the SD and the FD, are used to inter-calibrate 
these two detectors, providing an energy estimation which 
is almost independent of Monte Carlo simulations. Firstly, 
we have identified a SD observable which could be used to 
estimate the shower energy, and this was taken be the estimated 
signal 1000 m (S'(IOOO)) from the reconstructed shower core. 
Secondly, as such a signal is attenuated by the atmosphere, 
a correction is applied using the so called Constant Intensity 
Cut (CIC) [17]. After the correction, S'(IOOO) is normalized to 
a reference angle (38 degrees) and the final observable (6'38o) 
can then be correlated to the calorimetric energy measurement 
performed by the FD [18]. Once such a calibration curve is 
built for the hybrid events, it can be used for the whole high 
statistics sample measured by the SD. 



^Doublets are pairs of tanks deployed very close to each other (11 m), so 
that they sample essentially the same part of the shower. 



III. Establishing the Anisotropy of the Highest 
Energy Sample 

A. The source catalog and the scan method 

The Veron-Cetty (V-C) catalog [19] is a compilation of 
different surveys and its 12th edition contains 85221 quasars, 
1122 BLLacs and 21737 AGNs (of which 694 are less 
than approximately 100 Mpc from us). By considering only 
its closest sources {z < 0.024), the probability for chance 
correlations with sources from the catalog, already taken into 
account the Auger exposure, is as high as approximately 70% 
for angular windows of 7°. In other words, it is still possible 
to look for anisotropics. As farther sources are aggregated, of 
course, chance correlation becomes essentially 100%, making 
searches for deviations from isotropy impossible. 

Anisotropy signals in the cosmic ray arrival directions might 
be identified with the help of an intermediate source catalog 
through which a correlation signal might be established using 
an appropriate set of parameters such as: the typical angular 
scale of the correlation signal, the maximum redshift of the 
candidate sources and the minimum energy of the events. A 
scan over these parameters can then be performed looking for 
the typical values leading to a minimum of the probability for 
such a data-catalog correlation has been produced by chance. 
Given N events above a certain energy threshold Emin, and 
a set of sources (limited by a maximum redshift Zmax), the 
probability that at least k of them are at less than a certain 
angular distance 7 from one of the catalog's sources, just by 
chance, is given by the cumulative binomial probability 

^ / AT \ 

where pij^Zmax) is the probability for an event, drawn 
uniformly on the sky (but taking into account the detector 
exposure), to be at an angular distance less than 7 from at 
least one of the sources with z < Zmax- 

B. The exploratory scan 

A potential anisotropy signal in a subset of the Auger SD 
data, collected between 1 January 2004 and 27 May 2006, 
was identified through a scan over the 3 parameters mentioned 
above. Angular windows from 1° to 6°, V-C sources up 
io z = 0.024 and showers with energies above 40 EeV 
were considered. The minimum and maximum values for the 
windows are limited by the detector angular resolution and the 
constraint of moderate chance correlation probabilities with 
sources up to a certain maximum redshift, respectively. The 
soil energy was chosen in order to cope with the fact that 
typical magnitudes for the galactic and extra-galactic magnetic 
fields, will imprint deviations on low energy primaries, even 
protons, larger (in average) than the angular interval consid- 
ered. The minimizing parameters for this exploratory scan 
were: 7 = 3.1°, Zmax = 0.018 and Emin = 56 EeV. For 
such a configuration, 12 out of 15 events were correlated with 
V-C sources when only 3.2 were expected by chance (21% of 
chance correlation probability). 



C. The signal confirmation 

The Pierre Auger Collaboration has an internal policy which 
states that any potential anisotropy signal should be tested 
on an independent data set [20]. Therefore, a prescription 
was established in order to confirm or to reject the potential 
signal. Given the very low statistics at the energies considered, 
the Collaboration decided to proceed through a sequential 
analysis with a predefined stopping rule, since such tests can 
give an answer at earlier stages in the sequence. We have 
chosen isotropy as our null hypothesis, adopting a a = 1% 
Type I error, that is, the allowed probability for a false claim 
of anisotropy, and a Type II error of P = 5%, that is, the 
probability for an incorrect claim that the data was isotropic. 
The error (3 determines the power of the test (1 — (3), that 
is, 95%. By moving the bounds of the exploratory scan, the 
correlation power could drop from 80% (12/15) to values as 
low as 60% (8/14). In order to keep the test power above the 
95% threshold, even in the case of such a lower correlation 
power, the total length of the test was set to 34 additional 
detected events above 40 EeV (the predefined stopping rule). 

This prescription was fulfilled almost a year later on 25 
May 2007, with the detection of 8 new events above 40 EeV, of 
which 6 were at less than 3.1° from AGNs up to 75 Mpc from 
us. We could then claim that the hypothesis of isotropy for the 
highest energy Auger events (E > 40 EeV) was rejected with 
at least 99% CL. 

Hypothesis tests done via a classical sequential analysis 
have to deal with the fact that a penalty factor should be 
included at each intermediate step where the hypothesis has 
been tested, simply due to the fact that by performing the 
test several times, the probability for fulfilling it by chance is 
enhanced. Such a penalty factor will therefore depend on the 
total length of the test, introducing an undesirable dependence 
of the corresponding probability being monitored on data 
actually not taken. These classical tests drawbacks can be 
avoided by the use of a likelihood ratio test [21]. In figure 
[3] one can see the evolution of the likelihood ratio [22] 

as a function of time and the rejection of the isotropy hy- 
pothesis after the collection of 10 events out of which 7 were 
correlated with AGNs for the prescription parameters. 

D. Combining the full data set 

Once the anisotropy signal has been confirmed with the 
independent data set, all the data collected since 1 January 
2004 was submitted to the same scan, using as well a new 
calibration available at the time. The minimum was achieved 
for the 27 events above 57 EeV, of which 20 were at less than 
3.2° from AGNs of the V-C catalog up to approximately 71 

^The numerator is the probability to get a particular data set in case 
the anisotropy hypothesis is true and the denominator the corresponding 
probability when the isotropy hypothesis is true. Notice that since we do not 
know the true correlation power pi in case of anisotropy, we have integrated 
over this variable. 



Fig. 4. Sky map in galactic coordinates showing the directions of 27 Auger highest energy events (E > 57 EeV) with angular windows of 3.2 ° around 
them, as well as the positions of the 442 AGNs in the V-C catalog with redshift z < 0.017 (~ 71 Mpc). The color code indicates the observatory relative 
exposure over the sky. Centaurus A is marked in white. 
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Fig. 3. Likelihood ratio R (see eq.[2| evolution for every new event detected 
since the start of the internal Auger prescription. The threshold for isotropy 
rejection is represented by the dashed horizontal line. The 68% and 95% CL 
regions from isotropic expectations are also shown. Isotropy rejection was 
achieved after 10 detected events. 



Mpc (Zmax =0.017). Only 5.6 events are expected to correlate 
by chance under such conditions, and the probability for that 
to have happened by chance (already penalizing for the scan 
procedure) is approximately 10~^. 

Figure |4] contains a sky map in galactic coordinates with 
the 27 events above 57 EeV for the full data set with circles 
of 3.2° around each of them, as well as the positions of the 
442 AGNs with z < 0.017. The color code indicates the 
Auger relative exposure in each direction of the sky. A clear 
property of the events distribution is the concentration close 
to the supergalactic plane (dashed line), following the local 
anisotropic matter distribution in our neighborhood. Notice as 
well the two events at less than 3° from Centaurus A, one of 
the closest AGNs. 

IV. Signal Properties 

A closer look at the signal behavior in the vicinity of 
the minimum in the parameter space, shows the presence of 
additional secondary minima, indicating that the minimizing 
parameters shown above should not be taken at their face 
values, but only as indicative of the true correlation scales. 
Another distinctive property of the signal evolution as a func- 
tion of energy is the sharp transition to the global minimum. 



suggesting that something abrupt happens around E^r 
have the effect of enhancing the anisotropy signal. 



and 



Figure |5(a)| shows the average log-likelihood per event LL 
as a function of the energy threshold, defined as 



1 ^ 



(3) 



k=i 



where the sum runs over all the events above the energy 
threshold and p/c is the column density in the direction of 
the /c-event of a map built by smoothing out the V-C catalog 
with Gaussians of 2° around each of its sources. The plot 
shows also the la expectation bounds of the log-likelihood 
for sets of events (consistent with the statistics in the data 
for the corresponding energy threshold) drawn from the V- 
C catalog itself (red lines) and from an isotropic distribution 
(blue lines). One can clearly see the abrupt transition from 
isotropy to V-C-like anisotropy as the energy reaches 60 EeV 
from below. Such a sharp transition is accompanied by a 50% 
drop in the flux, from what would be expected from a simple 
extrapolation of the measured flux just below the cu toff, i n the 
form of 3. E~^-^^ power law as can be seen in figure 5(b) This 
striking coincidence is an additional evidence for the existence 
of the GZK effect, in which the sudden horizon shrinking, as 
the energy threshold for pion photo-production is reached, has 
the immediate consequence of revealing the anisotropics of our 
clumpy neighborhood. If sources were simply running out of 
power at energies around 60 EeV, such a sudden appearance 
of anisotropy would not be expected. 

Some idea of the anisotropic nature of the 27 most en- 
ergetic events detected by the Auger Observatory might be 
obtained without using the V-C catalog, therefore avoiding the 
incompleteness issues related to this compilation of sources. 
One of the most simple intrinsic properties of a given set of 
arrival directions is given by its auto-correlation function, that 
is, the number of pairs for different angular scales. The auto- 
correlation for the highest events identified in the scan over the 
full data set is shown in figure [6j together with the expectations 
from isotropic data sets. The departure from isotropy can be 
seen in a wide range of angles from 9° to 22°, where the data 
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Fig. 5. (a) Average log-likelihood per event (dots) for the Auger events as function of the energy threshold in EeV compared to the expectations from 
samples drawn from the V-C catalog smoothed at a 2° scale (dashed red line) and from isotropic samples (dashed blue line). The solid lines represent la 
bounds, (b) top: Differential flux as a function of energy and its statistical uncertainty. The number of detected showers for each data point is shown on top of 
each point; bottom: Fractional difference between the Auger cosmic ray SD spectrum and an ^-2-69 power law. Notice the coincidence between the onset 
of anisotropy for the events above ~ 60 EeV in the plot on the left-hand side and the energy where the flux drops by a factor 2 from what would be expected 
from a power law on the right-hand side. The HiRes data is also shown for comparison. 
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Fig. 6. Auto-correlation of the 27 highest energy events (dots) as a function 
of the angular scale in degrees. The expectations from isotropic samples of 
27 directions is also shown (crosses) and the error bars represent the 90% 
CL. 



shows a larger number of pairs than the isotropic expectations. 



V. Additional Remarks 

The topology of the galactic magnetic field is highly un- 
known. The few measurements of polarized light from galactic 
and extragalactic pulsars seems to point to a field of some /iG 
with direction reversals between arms and an exponentially de- 
creasing magnitude as a function of the distance to the galaxy 
center measured in the galactic plane [23]. With the increasing 
statistics of the Auger surface detector, the prospects for the 



detection of multiplets are very promising, which could bring 
valuable information on the galactic field structure. However, 
a true knowledge of the field topology should be provided by 
a dedicated experiment, such as the Square Kilometer Array 
(SKA) [24], which plans to measure about a million pulsar 
rotation curves, building a true 3 -dimensional map of the 
regular magnetic field component in our galaxy. This kind of 
information, combined with the detection of multiplets, would 
enable us to measure particles rigidity dependence with energy 
and, in turn, get information on their chemical composition. 

As already mentioned, the coincidence in energy between 
the flux attenuation and the enhancement of the anisotropy 
signal is an additional evidence for the GZK effect. Therefore, 
one could naively identify the maximum distance D^ax of 
sources coming from the scan with the horizon at ^min 
(defined as the distance from us containing the sources respon- 
sible for 90% of the protons above a certain energy threshold.). 
At a first sight, then, the value of 75 Mpc for Dmax as 
compared to ~ 200 Mpc [25] for the proton horizon at 60 EeV 
(assuming continuous energy losses, a uniform distribution of 
sources with equal intrinsec luminosities) could be seen as 
an inconsistency. Nonetheless, D^aax and the horizon are not 
even linearly correlated, chance correlations with foreground 
sources induce some bias on Dmax towards smaller source 
distances, and the same effect can be obtained in the presence 
of a large local overdensity. 

The HiRes Collaboration [26] has performed a similar 
analysis, but it does not seem to confirm the Auger result. 
Above 56 EeV, from a total of 13 events, only 2 are correlated 
with AGNs with z < 0.018. A difference in energy scale 



larger than the one envisaged by the Hires Collaboration 
could explain the discrepance if one believes the Northern 
and Southern skies are not so much different. On the other 
hand, an apparent lack of events in the region corresponding 
to the Virgo cluster seems to be a property of both data sets, an 
effect which could arise in a scenario where the bulk of the flux 
reaching us is generated in a very limited number of sources, 
like has been suggested in [27], in particular, analyzing the 
possibility of being Cen A the sole source. The construction 
of the Auger North [28] will allow us to crosscheck the 
anisotropy results obtained from the South Observatory and 
from HiRes, to enhance the statistics in the highest energy 
part of the spectrum, look for multiplets, and even measure 
individual source spectra. 

A detailed study of the sources correlating with the Auger 
high energy events, shows that they do not constitute the sub- 
sample with the most powerful sources. It is worth stressing 
though, that we do not claim the identification of the cosmic 
rays sources, and given the size of the correlation scale of less 
than ~ 6°, it is impossible to pinpoint the true ones, and es- 
sentially, anything whose spatial distribution, when convoluted 
with the magnetic deflections, follows the AGNs clustering 
properties, cannot be excluded so far. An interesting possibility 
raised in [29] is that, perhaps, some of the Auger events are 
pointing back to the last scattering centers, instead of their 
original sources, due to the presence of an inhomogeneously 
magnetized intergalactic medium. 

There have been some concerns with respect to the incom- 
pleteness of the V-C catalog, close to the galactic plane (due 
to dust extinction), at high redshifts and the lack of a selection 
function, since this is in fact a compilation of surveys and not a 
survey by itself. However, the demonstration of the anisotropy 
is not affected by these problems, since the catalog has been 
used only as an intermediate step in order to establish the 
typical scale of the anisotropy signal. 

VI. Conclusion 

The Pierre Auger Collaboration has established with more 
than 99% CL the anisotropy of its highest energy sample. 
The reconstructed directions of showers with energy above 
60 EeV show a correlation on angular scales of less than 
6° with the positions of AGNs up to ~100 Mpc from us. 
The evolution of the signal with the energy shows a sharp 
transition from isotropy to anisotropy at E ^ 60 EeV, and 
such a transition coincides with the flux attenuation observed 
in the Auger spectrum. This is an additional evidence for 
the existence of the GZK cutoff, where the abrupt reduction 
in the particles mean free path as the threshold for pion 
photo-production in the CMB is achieved, makes the universe 
essentially opaque at scales much larger than 100 Mpc for 
particles above 60 EeV. As a byproduct, the small effective 
horizon at high energies enhances the anisotropics signal, 
revealing the clumpy aspects of the structures in our local 
universe. Even though the correlation with the V-C catalog 
seems to be quite robust, the angular scale of ~ 6° does 
not make possible to unambiguously identify the sources. 



Essentially, anything which clusters (after taking into account 
magnetic deflections) in a similar way as AGNs are cannot be 
excluded as the true sources. Intrinsic (catalog independent) 
properties of the events, such as their auto-correlation function, 
show a clear departure from isotropy in a large angular range. 
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